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We present a comprehensive study of the evolution of the nematic electronic structure of FeSe using high 
resolution angle-resolved photoemission spectroscopy (ARPES), quantum oscillations in the normal state and 
elastoresistance measurements. Our high resolution ARPES allows us to track the Eermi surface deformation 
from four-fold to two-fold symmetry across the structural transition at ~ 87 K which is stabilized as a result 
of the dramatic splitting of bands associated with dxz and dyz character in the presence of strong electronic 
interactions. The low temperature Eermi surface is that of a compensated metal consisting of one hole and two 
electron bands and is fully determined by combining the knowledge from ARPES and quantum oscillations. 

A manifestation of the nematic state is the significant increase in the nematic susceptibility approaching the 
structural transition that we detect from our elastoresistance measurements on EeSe. The dramatic changes in 
electronic structure cannot be explained by the small lattice effects and, in the absence of magnetic fluctuations 
above the structural transition, point clearly towards an electronically driven transition in EeSe, stabilized by 
orbital-charge ordering. 


A nematic state is a form of electronic order which breaks 
the rotational symmetries without changing the translational 
symmetry of the lattice, and this state may play an impor¬ 
tant role in understanding high temperature superconductivity. 
There has been great interest in determining such electron ne¬ 
matic states m, most recently in iron-based superconductors 
la. Several scenarios have been proposed to explain the ob¬ 
served nematicity based on phonon-driven and electronically- 
driven tetragonal symmetry breaking. As the lattice effects 
are very small and such a state is found in close proximity 
to a magnetic state, it has been suggested that the nematic¬ 
ity is generated by spin fluctuations and is a precursor of the 
incipient antiferromagnetic state El. Other models empha¬ 
sizes the tendency towards orbital order for the Fe d elec¬ 
trons in different local environments, resulting in the devel¬ 
opment of unequal occupancy of dxz and dyz orbitals which 
drives the symmetry-breaking, and the magnetism is stabi¬ 
lized as a consequence m Electronic nematic order can be 
revealed through its strongly anisotropic resistivity and find¬ 
ings so far in iron-based superconductors have supported the 
spin-fluctuations scenario min. Understanding the dominant 
electronic interactions is of great importance, for they could 
determine the symmetry and properties of the superconduct¬ 
ing state 0. 

With a superconducting transition of Tc ^ 9 K, FeSe is a 
special case amongst Fe-based superconductors, since it un¬ 
dergoes a structural transition at ^ 87 K but does not order 
magnetically at any temperature. It has also attracted a lot of 
interest due to the strong increase in Tc to 37 K under pressure 
161 . the existence of high-Tc intercalates of FeSe (71 and, that 
a monolayer of FeSe grown on SrTiOs could have its super¬ 
conducting transition temperature in excess of 100 K (81. The 
availability of high quality bulk crystals, grown using chem¬ 


ical vapour transport El, have recently re-opened investiga¬ 
tions into the electronic properties of FeSe. ARPES has found 
evidence of a large band splitting caused by orbital ordering 
below the structural transition mini but the resolution of 
the available data cannot clarify the changes that occur at the 
Fermi level. Moreover, quantum oscillation experiments at 
low temperatures have detected an unusually small Fermi sur¬ 
face dlldSl. As magnetic fluctuations are detected only be¬ 
low the structural transition in FeSe, it is expected that they 
are not the driving force for this transition (TH [151. Thus, 
the nature of the structural transition in FeSe is rather unusual 
and how the electronic structure is stabilized by breaking of 
the rotational symmetry could be the key to understanding its 
superconductivity and how it can be further enhanced. 

In this paper we explore the electronic structure of FeSe by 
using high-resolution ARPES, quantum oscillations and elas- 
toresistivity measurements on the same batch of high-quality 
FeSe single crystals. These techniques provide a compre¬ 
hensive picture of the evolution of the electronic structure 
of FeSe from the high-temperature tetragonal phase through 
the fourfold-symmetry-breaking structural transition at Tg and 
into the electronic nematic phase. We observe strong in-plane 
(i-wave like deformation of the Fermi surface from four-fold 
to two-fold symmetry as a result of band splitting with dxz and 
dyz character. A direct manifestation of a nematic Fermi sur¬ 
face of FeSe, is the significant increase in the nematic suscep¬ 
tibility measured by elastoresistance measurements when ap¬ 
proaching Tg, not emerging due to anisotropic magnetic fluc¬ 
tuations but being driven by the orbital/charge degrees of free¬ 
dom. The low temperature Fermi surface, based on our high 
resolution data from ARPES and quantum oscillations, con¬ 
sist of an in-plane distorted quasi-two dimensional hole band 
and two electron bands; the inner electron pocket is extremely 
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elongated and quasi-two dimensional, whereas the outer elec¬ 
tron band, with predominantly d^y character, is not detected 
by ARPES but is present in quantum oscillations. The high- 
temperature bands have orbital dependent band renormalisa¬ 
tion, with respect to band structure calculations, due to many- 
body interactions, which are particularly significant for the 
dxy band; this leads to a significant shrinkage of the Fermi 
surfaces (a factor of ^5, compared with calculations). Our 
measurements also detect the band splitting induced by spin- 
orbit coupling 20 meV at the V point). 

Experimental details Samples were grown by the 
KCI/AICI 3 chemical vapour transport method. Magnetotrans¬ 
port measurements were performed as function of in-situ rota¬ 
tion in high magnetic field up to 33 T at HFML in Nijmegen, 
using an excitation current of 0.8 mA. Good electrical con¬ 
tacts were achieved by using In solder. Quantum oscilla¬ 
tions were observed in more than three samples with good 
agreement. ARPES measurements were performed at the 105 
beamline of Diamond Eight Source, UK. Single crystal sam¬ 
ples were cleaved in-situ in a vacuum lower than 2 x 10 “^^ 
mbar and measured at temperatures ranging from 6-120 K. 
Measurements were performed using linearly-polarised syn¬ 
chrotron light from 20-120 eV and employing Scienta R4000 
hemispherical electron energy analyser with an angular res¬ 
olution of 0.2-0.5 deg and an energy resolution of 3 meV. 
Elastoresistance measurements were performed by measuring 
the in-plane anisotropic transport properties while straining 
the crystals along the [ 110 ] direction in the tetragonal phase, 
by applying a voltage to a piezoelectric stack, similar to those 
discussed in Ref. 0 and detailed in SM. Band structure 
calculations were performed using the experimental crystal 
structure of FeSe in Wien2k using the GGA approximation 
with spin-orbit coupling included. Since FeSe does not show 
magnetic order, ARPES spectra are compared to non-spin- 
polarized band structure, with the relaxed lattice parameters, 
a = 3.7651 A, c = 5.5178 A, zse = 0.24128 CS). 

Temperature dependence of the hole bands. Fig. 
shows the calculated Fermi surface of FeSe in the tetrago¬ 
nal paramagnetic phase, which consists of three hole pockets 
{a, [3 and 7 ) around the Brillouin zone center (F) and two 
electron-like pockets ((5 and e), similar to other reports (Vfl . 
The bands which cross the Fermi level have mainly dxz^ dyz 
and dxy character (Fig[^ and b). Due to the matrix elements 
effects in ARPES experiments the light polarization and the 
scattering geometry allows us to select mainly either dxz or 
dyz bands using EH (p) or LV (s)-polarization, respectively, 
and to identify the orbital character of the measured bands, 
as shown in Supplementary Material (SM) and detailed in 
Ref.IISl 

Fig. [^-f show the evolution of the hole bands by perform¬ 
ing high-symmetry cuts, either through M-T-M {kz = 0) or 
A-Z-A (kz = tt / c ), as shown in Fig.[^-f as a function of tem¬ 
perature. In the high-temperature tetragonal phase with pre¬ 
served C 4 rotational symmetry, both a and p bands cross the 
Fermi level at the Z point, as shown by the Fermi-level mo¬ 
mentum distribution curves (MDGs) at 120 K in Fig[TJ. The 
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FIG. 1. Temperature dependence of the hole bands in FeSe. a) 

and b) Band structure calculation of the Fermi surfaee of FeSe in 
the tetragonal phase, projected into the kz = 0 plane and colored 
by dominant orbital charaeter; the high symmetry cuts are indicated. 
Calculations predict three hole poekets around the F point and two 
electron pockets around the M point, c) The high-symmetry M-F- 
M cut (cut 1) at low temperature (10 K). d)-f) The temperature- 
dependence of the A-Z-A cut at the Z point (kz = tv I c plane), g) and 
h) Extraeted band dispersion from constrained multiple Lorentzian 
fits to momentum-distribution eurves (MDC, solid symbols) or peaks 
in the energy dispersive cut (EDC, open symbols) from the data in d) 
and f). The solid lines in g) are renormalized band structure from 
a). The insets show a schematic Fermi surface, i) MDC for the hole 
bands (a and /3) at high and low temperatures. At 10 K, the a splits 
near the Fermi level, resulting in two crossings due to formation of 
twin domains, j) Fermi surface map at 10 K indieating the elongation 
of the hole pocket eombined with a duplicate rotated by 90° caused 
by twin domains, k) Photon-energy dependence of the MDC ai Ep 
(equivalent to scanning kz in the F-Z direction) that shows a quasi- 
two dimensional hole band at 10 K. 
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degeneracy between the a and P bands is lifted by the spin- 
orbit coupling insi which we can directly estimate from our 
data as being A^o ^ 20 meV, similar to results reported for 
LiFeAs 1^ . On lowering the temperature below the struc¬ 
tural transition at ^ 87 K (Fig[^, the [3 band at Z is 
pushed below the Fermi level, as shown in Fig. About 
50 meV below the Fermi level we also detect the presence 
of a broad band that is the dxy band, 7, in Fig. [^, having 
weaker intensity due to matrix elements and being shifted 
significantly in relation to the band structure calculations for 
the tetragonal structure in Fig. [^. The effect of the loss of 
C 4 rotational symmetry through the structural transition from 
tetragonal to orthorhombic symmetry at Tg is clearly observed 
in the detailed Fermi surface map at the Z point, showing that 
the a pocket elongates to become elliptical, shown in Fig[^ 
and j. Due to the natural twinning of the samples, however, 
we observe two ellipses (in Fig[T]i) rotated 90° with respect 
to each other. The kz dependence can be determined by mea¬ 
suring the Fermi-level MDC as a function of incident photon 
energy, shown in Fig. [^; we detect a quasi-two dimensional 
shape for the a hole pocket, the only hole band that crosses 
the Fermi level at low temperature. The experimental area for 
the hole bands, at 120 K is strongly reduced by a factor 
5, as compared with band structure calculations in the tetrag¬ 
onal phase (see SM) and we will later discuss the consistency 
with our quantum oscillations results at low temperatures (see 
Table 1 in SM). 

The effect of electronic correlations can be estimated by 
comparing band structure calculations with the band disper¬ 
sion in the tetragonal phase at 120 K, as shown in Fig. [T^. 
The obtained band renormalisation factors for FeSe are ^ 3.2, 
2.1 and 8 for the a,[3 and 7 pockets respectively - suggest¬ 
ing particularly strong orbital-dependent electronic correla¬ 
tions on the 7 band with d^y orbital character. These values 
are close to those reported for FeSe in Ref. ini, whereas in 
the case of FeSe^^Tei-a^, the band-selective renormalization 
varies strongly with values between 1 and 17 1^ . the sig¬ 
nificantly higher value being assigned to the dxy band. Band 
structure calculations including correlations (DMFT) on FeSe 
II 22 I give band renormalisations of ^ 2.8 for the dxzjyz or¬ 
bitals, comparable to the measured values, whereas the pre¬ 
dicted value of ^3.5 for the dxy band is smaller than in exper¬ 
iments. 

Temperature dependence of the electron-like bands. We 

now focus on the temperature dependence of the electron¬ 
like bands in FeSe shown in Fig. [^, visualizing the dramatic 
changes occurring through the structural transition. In the 
high-temperature tetragonal phase of FeSe, band structure cal¬ 
culations predict two cylindrical electron pockets, as shown 
in Fig.[^ and Fig.[^. The outer electron-like band (labelled 

has a dominant dxy character whereas the inner electron 
band (labelled e) has dxzjyz orbital character. At high temper¬ 
atures (120 K), the inner dxzjyz band has a strong intensity 
at the M point, where bands disperse both along and perpen¬ 
dicular to the cut direction and are four-fold degenerate, as 
shown in Fig|^. However, a combination of low intensity 


caused by matrix element effects for the dxy orbital (HI, as 
well as broadening due to either impurity scattering 12 ^ or 
strong correlations may cause band incoherence (221 1231 and 
difficulty to observe the dxy band in ARPES at the M point 
CD ESI [261. 

Following the temperature dependence of the inner elec¬ 
tron bands with dxzjyz orbital character, in the experimen¬ 
tal F-M-F cut (cut 2 in Fig[^), the band crossing the Fermi 
level would have dyz character, whereas the downwards dis¬ 
persing band from the M point would have dxz character (see 
Fig.|^). Below the structural transition which breaks the C 4 
rotational symmetry, we observe a large band splitting which 
can only occur if the dxzjyz degeneracy is broken and the dxz 
and dyz orbitals develop an unequal occupation. As a result, 
the dyz band moves up (ei in Fig.[^), and the dxz band moves 
down (solid lines), but for the other structural domain the op¬ 
posite occurs, the dxz band moves up and the dyz ( 62 ) moves 
down indicated by the dashed lines in Fig. and further de¬ 
tailed in SM {x and y are defined in the experimental coor¬ 
dinate frame). This picture is consistent with that suggested 
from ARPES studies on detwinned crystals of FeSe (TH . but 
the high resolution of our data allows us to detect the pre¬ 
cise changes in the Fermi surface across the transition. The 
large band splitting of ^ 50 meV at the M point at 10 K 
(Fig. 1^) indicates the lifting of of dxz and dyz degeneracy 
in FeSe. This was also observed in NaFeAs and BaFe 2 As 2 
over a very limited temperature range between the structural 
transition and magnetic transition, |[T^[27H29ll . In FeSe, this 
dramatic energetic shift observed below Tg is much larger than 
that expected from a simple structural orthorhombic distortion 
of 2 X 10“^ ((a — h)l{a + h)), which would cause a 17 meV 
shift in the absence of the renormalization effects (see SM). 
At the Fermi level, the resulting electron pocket is strongly 
elongated, with a kpiv) ~ 0.02 and kpix) « 0.14 
A“^. However, due to twinning, the Fermi surface consists of 
two ellipses at 90° to each other, resulting in a cross-shaped 
Fermi surface (Fig. |^, similar to the hole pocket. The de¬ 
gree of Fermi surface distortion is measured by the tempera¬ 
ture dependence of kp, and the ellipticity kF{y){si), shown 
in Fig. 1^, is suggestive of an order parameter of a second or¬ 
der phase transition at Tg. Thus, the observed experimental 
elongation of both the electron and hole pockets at low tem¬ 
peratures in FeSe can be thought of as a consequence of the 
electronic anisotropy induced by orbital ordering in the pres¬ 
ence of interactions. Another possible scenario for this be¬ 
haviour can be thought as a Pomeranchuk instability (301 that 
results in the spontaneous deformation of the Fermi surface 
from a four-fold symmetric almost circular shape at high tem¬ 
perature to an elliptical shape at low temperatures. In either 
cases, this electronically-induced Fermi surface deformation 
is also supported by the small specific heat jump at Tg com¬ 
patible with an electronic contribution (9l. 

Comparison with quantum oscillations. In order to pro¬ 
vide a complete and consistent picture of the Fermi surface of 
FeSe, we have also measured quantum oscillations in samples 
from the same batch used in the ARPES experiments. Quan- 
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FIG. 2. Development of a strongly elongated electron Fermi surface through the structural transition, a) Temperature dependence of 
high-symmetry cut through M, showing the development of a splitting of intensity which results from strong orbital-dependent shifts of the 
band structure with onset at Ts = 87 K. b) The calculated band dispersions in the tetragonal structure, c) The experimental band structure 
at 120 K (solid dots) renormalized to the high temperature data (solid lines); the outer dxy electron band ((5) is not observed at all at the M 
point, d) The lifting of dxzjyz degeneracy at 10 K and the consequent band shifting for two different domains, ei and 62 , indicated by the 
solid and dashed lines, respectively. The insets show the schematic band dispersion for different domains and directions, as detailed in SM. 
e) Temperature-dependence of the EDC at the M point, showing the development of the splitting of bands with dxz and dyz character; the 
splitting A is plotted in Fig. I^n. f) Fermi surface map (integrated within 2 meV of the Fermi level) showing the cross-shape arising from 
strongly elongated Fermi surfaces at M in the two domains, shown schematically in the inset. 


turn oscillations are a powerful technique that allows precise 
determination of cross-sectional areas of the Fermi surface 
and the corresponding orbitally-averaged quasiparticle masses 
but it usually needs to rely on band structure or ARPES to pro¬ 
vide the exact /c-space location of these orbits. Magnetotrans¬ 
port measurements are performed in the normal state of FeSe 
at very low temperatures and high magnetic fields, as shown 
in Fig. [^. The oscillatory signal, periodic in 1/B, is made 
clear by subtracting a high order polynomial from the raw 
data in Fig. [^. The fast Fourier transform of pose shown in 
Fig#, identifies four different quantum oscillation frequen¬ 
cies, labelled Fi _4 in ascending order, that are directly linked 
to the extremal areas on the Fermi surface by the Onsager re¬ 
lation, Fi = Akihl{27Te), for a particular field orientation. 
All frequencies are lower than 1 kT in agreement with previ¬ 
ous reports |[T2l [TSl . This is smaller than the frequencies of 
other iron-based superconductors that do not show any Fermi 
surface reconstruction, such as LaFePO ED or LiFeAsEa. 
Band structure calculations predict the existence of five dif¬ 
ferent quasi-two dimensional cylinders with sizes (1-2.6 kT), 
much larger than those found in experiments (below 1 kT) 
(see SM). The angular dependence of these frequencies also 
suggest a quasi-two dimensional nature of these bands (see 
that the angular dependence of frequencies varies almost like 
1 / cos 0 in Fig|^, and previous studies have suggested that at 
low temperatures FeSe has either only quasi-two dimensional 


electron bands or electron and hole bands 

The cyclotron-averaged effective masses of the quasiparti¬ 
cles, extracted from the temperature-dependence of the quan¬ 
tum oscillations amplitudes (Fig|^) for each orbit, are listed 
in Table 1 in SM). We find a reasonably good agreement be¬ 
tween the values of the effective masses for our different sam¬ 
ples with those reported previously in Ref.[T2]and[T3l The F 2 
and F 4 orbits have similar effective masses of around 4 irie, 
which may indicate that they originate from the same band, 
but with different -values, as in the case of a corrugated 
quasi-two dimensional band, whereas the F 3 orbit has a par¬ 
ticularly heavy effective mass of ^ 7(1) rUe, in agreement 
with other reports ifT^ . However, in the case of the Fi pocket 
there is a variation of the quasiparticle masses between 0.7(2) 
ca to 3(1) rrie. The strong disparity between the lighter 
masses Fi and the heavier masses of F 4 points towards dif¬ 
ferent band origin for these orbits at the Fermi level. 

Next, we compare the absolute size of the different k val¬ 
ues at the Fermi level extracted from ARPES with those from 
quantum oscillations (Table 1 in SM). At low temperatures, 
ARPES detects a single quasi-two dimensional hole band (see 
Fig. Ik) along the T(Z) direction, which has a two-fold sym¬ 
metry, with a maximum area around Z and minimum around 
r, and a small carrier density of 2.53 x 10^^ cm“^. Similarly, 
ARPES gives clear indication of the presence of an electron- 
band centered at M which is a strongly elongated cylinder, 
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FIG. 3. Quantum oscillations in FeSe. a) Magnetoresistance of 
FeSe for a field applied along the c axis (^=0) with the irreversibil¬ 
ity field, Birr around 14 T at 0.36 K. b) The oscillatory part of the 
resistivity A pose , obtained by subtracting a polynomial background 
from the raw data, c) Fourier transform (FFT) of the data in b) identi¬ 
fying four different frequencies Fi_ 4 , corresponding to /c-dependent 
extremal areas on the Fermi surface, d) Temperature-dependence of 
the amplitude of oscillation, from which the effective masses may be 
extracted, e) Quantum oscillations and f) the corresponding Fourier 
transforms as a function of angle 0 of applied magnetic field with 
respect to the c axis at 0.36 K. By plotting F cos 0 against 0 one can 
identify if orbits correspond to a maximum or minimum of a quasi- 
two dimensional Fermi surface. Color bar corresponds to the Fourier 
transform amplitude. 


but it is lacking the information concerning the behaviour 
of the outer electron-like d^y band crossing the Fermi level, 
which may be present according to band structure calcula¬ 
tions. Since ARPES data suggest that there is only one hole 
pocket but possibly more than one electron pocket, the largest 
maximum orbit observed in quantum oscillations, F 4 , must 
corresponds to the hole pocket in order to maintain charge 
balance in the system; thus two quantum oscillation frequen¬ 
cies with similar effective masses, F 2 and F 4 , belong to this 


quasi-two dimensional hole band, in close agreement with the 
size of this quasi-two dimensional band detected from ARPES 
(Eig[^). Eurthermore, the F 3 frequency in quantum oscilla¬ 
tions has a substantially larger effective mass than all the other 
frequencies, which is of similar value to the band renormal¬ 
ization found in the ARPES experiments for the dxy band (7 
in Eig[^), suggesting that the outer electron pocket that is de¬ 
tected in quantum oscillations but not in ARPES has dxy char¬ 
acter. The inner elongated electron band from ARPES (with 
areas smaller than 100 T), can be assigned to the Fi frequen¬ 
cies, as this band is almost two-dimensional (as in seen at the 
A-point in Eig.l in SM). Thus, by combining the knowledge 
about the sizes and quasiparticle masses from the ARPES and 
quantum oscillations data, we suggest that bulk EeSe at low 
temperatures has one hole and two-electron bands, as repre¬ 
sented schematically in Eig|^ There are some small discrep¬ 
ancies in the absolute values of the cross section areas ex¬ 
tracted from ARPES and quantum oscillations, but these could 
be caused by field-induced Eermi surface effects in quantum 
oscillations or possible surface effects in ARPES. This inter¬ 
pretation is also consistent with the three-band description of a 
recent magnetotransport study on EeSe 13^ . in which a very 
small high mobility carrier was also detected, in agreement 
with the size of small electron band from our ARPES data. It 
is important to note that our measured electronic structure of 
bulk EeSe is similar to that of multilayers of EeSe on SrTiOs, 
in which the splitting of bands at the M point was wrongly 
assigned to a possible formation of a spin-density wave 1261 . 
However, the EeSe monolayer is rather different and its elec¬ 
tronic and high Tc superconducting state was suggested to be 
influenced by the SrTiOs substrate 1^ . 

One important flnding from our the experimental studies 
of the electronic structure of EeSe is that the band struc¬ 
ture calculations signiflcantly overestimate the size of the 
Eermi surface of EeSe, even in the tetragonal phase. Shrink¬ 
ing of Eermi surfaces in other iron-based superconductors 
has been assigned to the interband coupling to a bosonic 
mode in LaEePO, Ea, and/or the strength of the antifer¬ 
romagnetic correlations close to a quantum critical point in 
BaEe 2 (Asi_a;Pa ^)2 1^ . It is clear that in EeSe, the signifi¬ 
cant shrinking of the Eermi surface, combined with the strong 
renormalization effects both in ARPES and quantum oscilla¬ 
tions as evidenced by the relatively large effective masses of 
4-8 rrie and the orbital-dependent correlations (largest for the 
dxy band) suggest that the electronic correlations signiflcantly 
affect the electronic structure of EeSe. 

The nematic susceptibility of FeSe. A clear manifestation 
of a nematic Eermi surface is its strong in-plane anisotropy in 
transport properties and sensitivity to external parameters, in 
particular in-plane strain. The resistivity anisotropy is deter¬ 
mined by both the electronic structure and the scattering, and 
the expected Eermi surface deformation give rise naturally to 
anisotropic electronic properties, whereas the spin-nematic or¬ 
dering leads to an anisotropy of the electron scattering 0 . 

Eig|^ shows the induced change in resistivity to in-plane 
strain measured by the 2 m 66 parameter from elastoresis- 
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FIG. 4. Symmetry-breaking of the electronic structure of FeSe. a) Band-structure calculation of the Fermi surface in the high-temperature 
PA/nmm tetragonal space-group, b) Schematic and c) the experimental high-temperature in-plane Fermi surface consists of significantly 
shrunk pockets, d) Elongated hole and electron pockets below the structural transition, e) Schematic and f) experimental low temperature 
Fermi surface including the effect of twin domains (the dxy electron pocket is not plotted here), g-i) Three-dimensional representations of 
Fermi surface of FeSe, as described in a, b and d, respectively, j) Resistivity as a function of temperature for the sample used in Fig. 3. k) The 
induced change in resistivity to in-plane strain measured by the mee parameter from elastoresistance tensor in FeSe, which provides a direct 
measure of the electronic nematic order parameter, indicating that the structural transition is electronically-driven, as discussed in the main 
text. 1) Temperature-dependence of the intense kpi for the e pocket around the M point, showing the Fermi surface deformation which onsets 
at Ts. m) Energy splitting of the bands with dxz and dyz band character at the M point, extracted from Eig. 2e. 


tance tensor in FeSe, which provides a direct measure of the 
electronic nematic order parameter, as detailed in Fig.SMS 
and in Ref|5] We observe a large increase in 2m66 ap¬ 
proaching Ts, similar to large divergent behaviour observed 
in Ba(Fe/Co) 2 As 2 |[5][371, but of even larger magnitude. The 
data on FeSe can be well-described by a fit to the function 
2m66 = AjiT — T*) + Ao, which gives T* = 66(1) K and 
^qzz- 31(3) (solid line in Fig|^). Recently, nematic suscepti¬ 
bility measurements of the elastic shear modulus, which probe 
the lattice response to strain, suggest that the structural transi¬ 
tion in FeSe is accompanied by a large shear-modulus soften¬ 
ing, identical to that of underdoped Ba(Fe,Co) 2 As 2 , implying 
a very similar strength of the electron-lattice coupling, which 
is much smaller than the electronic response ifT^ . 

The sign of mee is also opposite to what is found in the 
electron-doped pnictides, but similar to FeTe m, where it 
was suggested that the resistance along the a (AFM) direction 
is larger than that along b axis (FM direction) (TSl. However, 


a small sign-change of the in-plane anisotropy was also found 
for highly hole doped (Ba/K)Fe 2 As 2 (Ml- This may suggest 
that the positive sign of mee may be either a general feature of 
chalcogenides, as opposed to most of the pnictides, or it is also 
possible that the development of the anisotropic properties in 
all these systems may be driven by a different mechanism. 

Below the structural transition, T^, the behaviour of mee 
is rather striking, having an almost linear temperature depen¬ 
dence with no sign of saturation, as the degree of elliptic- 
ity grows larger (see Fig|^); it changes sign around 65 K, 
which, coincidentally, is the same scale as the value of T* 
determined earlier. While the interpretation of nematic sus¬ 
ceptibility below Ts may be difficult due to domain formation 
and the fact that the nematic order parameter now takes a fi¬ 
nite value, we suggest that one possible explanation for this 
crossover may be linked to the development of anisotropic 
scattering from spin-fluctuation, which become strong below 
Ts cma. Sign-change of the in-plane anisotropy has been 
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found between the electron and hole-doped BaFe 2 As 2 ll38]| . 
being assigned to differences in the spin fluctuations scatter¬ 
ing rates corresponding to different Fermi velocities at the hot¬ 
spots for electron- and hole-doped systems |[38l[39l. Magne¬ 
totransport studies in FeSe also suggest that anisotropic scat¬ 
tering may develop below Tg 13^ . 

As the size of our measured nematic susceptibility is much 
larger than the response of the lattice (141, we suggest that the 
structural transition in FeSe is electronically driven and is an 
instability of the electronic structure which breaks tetragonal 
symmetry, with the lattice orthorhombicity simply responding 
to these electronic changes. Furthermore, the absence of spin- 
fluctuations above Tg indicates that the structural transition in 
FeSe is not magnetically-driven M\IB- Thus, our resistiv¬ 
ity anisotropy measurements favor an orbital/charge ordering 
scenario, which show a strong splitting of the bands with d^z 
and dyz orbital character, (Fig|^ and summarized in Fig|^) 
and it is likely responsible for the in-plane Fermi surface de¬ 
formation which give rise to signiflcant anisotropy in the in¬ 
plane electronic structure. 

While the orbital ordering scenario is a likely contender to 
explain the existing data on FeSe, the Fermi surface defor¬ 
mation and its relatively constant volume conservation (small 
changes occur due to the temperature dependence of the lattice 
parameters) bears the similarities of a d-wave Pomeranchuck 
instability. An isotropic Fermi liquid in the presence of suffi¬ 
ciently strong interactions was predicted by Pomeranchuk to 
be unstable fSOl and as a signature of this instability the Fermi 
surface would spontaneously deform, changing its shape or 
topology to lower its energy and break the rotational symme¬ 
try, similar to what is observed in FeSe at the Fermi level. The 
two scenarios are likely to generate different kinds of pairing 
interactions either at finite momentum ioi or at zero momen¬ 
tum ED and future theoretical work will need to address these 
issues. 

In summary, we report a comprehensive study of the de¬ 
velopment of the nematic phase in FeSe. The Fermi surface 
of FeSe undergoes a spontaneous distortion from fourfold- 
symmetric to two-fold symmetric elliptical pockets. The sym¬ 
metry breaking arises from the electronic degrees of freedom 
and, in the absence of magnetism, it is likely to be caused 
by orbital ordering in the presence of strong interactions. 
The elongated Fermi surface causes strongly anisotropic elec¬ 
tronic properties and enhanced nematic suceptibility. This ne¬ 
matic electronic phase is that from which superconductivity 
emerges, which is suggested to have a two-fold gap symmetry 
14^ . While interactions favoring orbital ordering dominate 
near Tg, magnetic fluctuations grow towards Tc, and may still 
assist the superconducting pairing. Moreover, the existence of 
a relatively flat hole band just below the Fermi level at the M 
point raises the question if whether the presence of van Hove 
singularities could also play a role in pairing. Given the ob¬ 
served small and strongly anisotropic Fermi surfaces with low 
carrier densities, it is perhaps not surprising that the physical 
properties, including the superconducting and structural tran¬ 
sitions are susceptible to external parameters (e.g. pressure. 


strain, doping) so that controlling these one could turn FeSe 
into a high-temperature superconductor. 
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